Resonances in dissipative optomechanics with nanoparticles: 
Sorting, speed rectification and transverse cooling 
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The interaction between dielectric particles and a laser-driven optical cavity gives rise to both 
conservative and dissipative dynamics, which can be used to levitate, trap and cool nanoparticles. 
We analytically and numerically study a two-mode setup in which the optical potentials along the 
cavity axis cancel, so that the resulting dynamics is almost purely dissipative. For appropriate 
detunings of the laser-drives, this dissipative optomechanical dynamics can be used to sort particles 
according to their size, to rectify their velocities and to enhance transverse cooling. 



I. INTRODUCTION 

In recent years, optically levitated dielectric 
nanospheres [U [2] have attracted considerable in- 
terest as a novel and versatile alternative to more 
conventional optomechanical setups [3HS]- A dielectric 
object inside a cavity with fixed mirrors gives rise to 
dispersive coupling between the center-of-mass motion 
of the object and the cavity mode, similar to the 
interaction between a cavity mode and a moving mirror 
[7] . Optically [IJ [2] or magnetically [8] levitated particles 
have the additional advantage that they are naturally 
well isolated from their environment. Combined with 
laser trapping and cooling, this may allow for cooling 
close to the quantum-mechanical ground-state [9], which 
offers the prospect to reach an entirely new regime in 
the ongoing pursuit to create quantum superpositions 
of particles with large masses [TO] E] and to study and 
probe unexplored regimes of quantum mechanics and 
possible deviations from it [HI H3] . 

Laser cooling in optomechanical systems has been 
studied extensively in the resolved-sideband regime (141 — 
[17] . which is analogous to the laser cooling of trapped 
ions [IB], and may be applied to cool molecules [T5]. In 
this regime, the nanoparticle is harmonically confined to 
half-a-wavelength of an optical trapping mode and laser- 
cooled through a second cavity mode pj [20] . Recently, 
the interaction between a nanoparticle and a cavity mode 
has been considered in a more general form [STjrjJ]. How- 
ever, laser cooling, velocity-dependent optical forces and 
dissipative dynamics in the regime in which the particle 
is not trapped but rather moves almost freely along the 
cavity axis, have not been explored yet. 

In this paper, we consider a setup with a nanosphere 
inside an optical cavity with two laser-driven modes, see 
Fig. [ija). The optical modes create a strong potential 
minimum along the cavity axis z. We consider two cases: 




FIG. 1: (a) A sub-wavelength nanoparticle of mass m is con- 
fined to the axis (z) of an optical cavity but can move freely 
along it. It is optomechanically coupled to two laser-driven 
cavity modes. The strengths are denoted by f2i,2 and the de- 
tuning from the cavity frequency is Ai,2. The decay of the 
cavity is k. (b) Effective damping force for one drive detuned 
to the red (see main text for the parameters used). For small 
velocities, the effective force exerted on the nanoparticle de- 
pends linearly on the velocity (dotted line) while resonances 
appear for larger velocities, (c) By detuning one drive to the 
red and the other to the blue with an additional relative de- 
tuning Ai = — A2 + S, the competing damping (dotted red) 
and heating (dotted blue) give rise to an effective dissipative 
force (black) with finite stationary velocities. Remarkably, 
since the force has a negative slope around these points, the 
stationary points are stable and all particles end up with the 
same final speed (determined by the detunings). 
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one in which the particles are injected in the cavity so 
that they will be trapped and move only along the cav- 
ity axis and one in which the particles are shot through 
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the cavity along the x or y axes and the velocity com- 
ponent along the cavity axis z is optically cooled during 
the passage time. In both cases, the spatial directions 
remain uncoupled and the dynamics is effectively one di- 
mensional. We choose the modes and drive strengths 
such that the conservative optical forces along the cav- 
ity axis almost completely cancel throughout the focal 
region of the cavity, and the resulting dynamics of the 
nanoparticle is dominated by dissipation. 

Laser cooling in this setup is reminiscent of Doppler 
cooling of moving atoms |23) . but it is also fundamen- 
tally different. An atom has discrete internal levels, 
which makes the effective force exerted by a laser beam 
strongly depend on the optical frequency and, through 
the Doppler effect, on its speed. Dielectric nanoparti- 
cles, on the other hand, do not have any resolved in- 
ternal structure and, in our setup, dissipation is purely 
due the finite cavity response time, which leads to a 
retarded, and, therefore, velocity-dependent, radiation- 
pressure feedback force. As one may expect, red-detuned 
laser drives (with respect to the bare cavity resonances) 
lead to damping (cooling), while blue detunings give rise 
to heating. 

In a two mode setup this naturally leads to two dif- 
ferent regimes: (i) in which the particle motion is cither 
damped or amplified by the optical modes and (ii) in 
which the motion is damped by one mode and amplified 
by the other, see Fig. [TJb)-(c). 

In case (i) with red detunings and in the limit of small 
velocities, the effective damping force depends linearly on 
the particle momentum, see Fig. [ljb), and gives rise to 
tunable spatially non- uniform Langevin dynamics |24j . 
For larger momenta, we find that the effective dissipa- 
tive force has a much richer structure with resonances 
that give rise to strongly enhanced damping, at least 
for specific momenta. The range of velocities for which 
such laser cooling works is significantly enhanced by the 
cancellation of the optical potential. We show that this 
can be relevant for quantum interference experiments in 
which nanoparticles are shot through a cavity and the 
transverse motion is optically cooled during the passage 
time [TO]. In this setup, transverse cooling is enhanced 
by almost two orders of magnitude. 

In case (ii), the dynamics is governed by the interplay 
between heating and cooling, see Fig. [TJc). If the two 
detunings are equal but opposite and the drive strengths 
are equal, the effects cancel and the motion is frictionless. 
However, if the magnitude of the detunings is (slightly) 
different, the non-linear nature of the effective forces lead 
to finite stationary velocities. For properly chosen de- 
tunings, these points correspond to stable stationary so- 
lutions. Initially slower particles are accelerated, while 
fast ones are decelerated to the same final speed. The 
effective dissipative force depends on the particle mass, 
but the stationary points in case (ii) do not. Using (ii) 
to initialize the particles at the same speed and, subse- 
quently, (i) to achieve mass-dependent damping, dielec- 
tric particles can be weighted and sorted according to 



their mass. In this setup, a mass difference on the order 
of Am ~ 10 -19 kg translates into a position difference 
on the order of Az ~ 10~ 3 m. 

This paper is organized as follows. In the next section, 
we introduce the setup in full detail and derive the equa- 
tions of motion. Numerical results for the full model are 
discussed in Sec. Ill, while an approximate analytical 
treatment is presented in Sec. IV. The analytical results 
are in excellent agreement with the numerical simula- 
tions and provide additional insights in the nature of the 
dissipative optical force and its possible application to 
velocity-rectification, sorting and laser cooling. In Sec. 
V, we discuss these application in more detail and a brief 
summary and some concluding remarks are contained in 
Sec. IV.' 



II. MODEL AND HAMILTONIAN 

We consider the setup shown in Fig. [lja), in which 
a sub-wavelength dielectric particle interacts with two 
laser-driven cavity modes. For moderately strong drives, 
the particle is confined to the cavity axis (z) and we 
assume that the variation of the mode functions in the 
transverse directions (x and y) is negligible compared to 
that along the z axis. In this case the dynamics is effec- 
tively one dimensional and, in a frame rotating with the 
drive frequencies, the classical Hamiltonian is given by 

U =t n - E U^la^-Q^+a*)) , (1) 

.7=1,2 V 1 

where j runs over the cavity modes, a\ t 2 are classical 
normal variables characterizing their dynamics, 0,1^2 are 
the drive strengths and Ai^(z) — ufl — wi :2 (z), with 

Lu[ d l the drive frequencies, are the detunings from the 
mode frequencies. The optomechanical interaction be- 
tween the particle and the cavity modes, derives from a 
dependence of the mode frequencies on the particle po- 
sition z along the cavity axis. We focus on the case of 
two modes separated by the a free spectral range of the 
cavity Au; = \lo^ — oj^\ = ctt/L, where L is the cavity 
length and uyf\ are the frequencies in the absence of the 
particle. We assume that Aw is much smaller than the 
absolute frequencies ui j°2 , but much larger than any other 
time scale involved. In this case the mode functions near 
the focal region of the cavity can be approximated by 
Fi(z) ~ e lv / 2/Csm(kz) and F 2 (z) ~ e 2 \/VC cos(kz) 
respectively. Here C — (-k/A)Lw 2 is the mode vol- 
ume, with w the beam width, and are polariza- 
tion vectors in the xy plane. The resonance frequencies 
are modulated by the intensity at the particle position 
(see appendix) so that oj\{z) — — gism 2 (kz) and 

uj\(z) = ^2°' — 52 cos 2 (fcz), where gj — uij°^a p V/C, with 
V the volume of the particle and a p its electric polar- 
izability, the optomechanical coupling strengths. Since 

Au> << cu'fl we assume that gi^ = g- It follows that 
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the position-dependent detunings are given by Ai(z) = 
Ai - (g/2)cos(2fcz) and A 2 (z) = A 2 + (g/2) cos(2/cz), 



(a) 



'1,2 



w i°2 T 5/2 are the average detun- 



where Ai )2 
ings. 

The decay of the optical modes is included at the level 
of the equation of motion hj — (iAj(z) — K,)cij — iQj. 
The drive term can be eliminated by moving to a shifted 
frame aj — » ctj + a'p where ctj — Clj/(Aj + in) is the 
average amplitude in the cavity modes. The equations 
of motion for the residual optical amplitudes a' are given 
by 

o-j = i iA j - k)clj + i(Aj(z) - Aj)(aj + a 3 ) , (2) 

where we have dropped the primes. The average mode 
amplitudes also enter the Hamilton equations for the par- 
ticle motion, which are found as 



P 



and 



2 dAj(z) 



(3) 



The \ctj\ 2 terms do not depend on the optical ampli- 
tudes dj and describe and optical potential V{z) — 
g{\oc\\ 2 — I ct 2 1 2 ) sin(2fcz). For g > 0, the particle ex- 
periences a force toward the high intensity regions of 
the optical field. This is the basic principle underlying 
optical tweezers [25]. The optical potential V(z) obvi- 
ously vanishes for equal intensities |ai| 2 and |a;2 1 2 • Since 
|(X,| 2 = |£X/| 2 /(A 2 + K 2 ), cancellation of the potentials 
requires that the drive strengths O12 depend on the de- 
tunings Ai j2 . Such cancellation is only exact for the 
approximate mode functions F 12 (^) introduced above. 
This is always an approximation and, in a typical setup, 
the strength of the effective optical potential depends on 
z according to cos 2 (ttz/ L), which vanishes near z ~ L/2. 
The potentials only cancel along the cavity axis. In the 
x and y directions, the optical fields create a potential 
minimum of — g(|ai| 2 + |o:2 1 2 ) along the cavity axis. The 
other terms in Eq. (|3| involve aj and couple the parti- 
cle motion to the dynamics of the optical modes. The 
cavity decay k sets a finite time scale for the cavity field 
to adopt to changes in the particle position, so that the 
radiation-pressure force on the particle is retarded and 
not only depends on the particle position, but also on 
the momentum along the cavity axis. Depending on the 
sign of the detuning, such forces give rise to damping 
(A < 0) or heating (A > 0). We focus on the case in 
which |ai| = |a 2 | so that, at least to leading order, the 
optical potentials cancel, and the dynamics is dominated 
by dissipation. 



III. NUMERICAL RESULTS 

Before we solve the equations of motion ([2| and ^ an- 
alytically by making further assumptions, we first study 
the full dynamics numerically. We assume that a single 
particle is carefully injected in the cavity so that it is 
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FIG. 2: (a) Velocity as a function of time for various 
nanosphere sizes ranging from r — 46 nm (black) to r — 54 nm 
(orange), (b) Position as a function of time for the same 
nanospheres. All particles are eventually trapped, but the 
distance traveled strongly depends on the size, (c) Zoom in 
to the point where a particle is trapped in the small reminis- 
cent optical potential. The result from the effective force Eq. 
(Tf| (black), the analytic result from Eq. Q (green) and the 
full numerical data (green) are in excellent agreement up to 
the point where the particle is trapped, (d) The amplitude of 
the cooling mode |ai|/|a| as function of time for R = 46 nm. 
The regime of strongly enhanced cooling in (a) corresponds 
to a sharp resonance in the optical response. 



trapped and, initially, moves along the cavity axis at a 
constant speed v(0). We numerically solve ([2| and Q for 
the z(t) and v(t) — p(t)/m for different particle masses. 
For concreteness, we consider Silica (p = 2200 kg/m 3 and 
e r = 2.34) nanospheres with radii around R = 50 nm, 
which have masses around m — 1.15 x 10 -18 kg. The 
optomechanical coupling strength g is proportional the 
particle volume (or mass) and, for a cavity length of L — 
0.12 m and a mode width of w = 100 /im, we find optome- 
chanical coupling strengths around g ~ 10 3 Hz. Further- 
more, we choose a wavelength A = 10~ 6 m for both modes 
and a cavity finesse of T = 50 x 10 3 , which corresponds 
to a field decay rate of k = ttc/(2LF) ~ 8 x 10 4 Hz. 
The intra-cavity power is taken to be P cav — 50 W 
per mode so that |ai. 2 | ~ 3.2 x 10~ 12 \fJs. For a 
50 nm particle, the optical confinement to the cavity axis 
— s(l a i| 2 + l a 2| 2 ) corresponds to a velocity of 0.19 m/s. 

For these parameters and with the first cavity mode 
detuned to the red Ai = —2k ~ —1.6 x 10 5 Hz and the 
other driven resonantly A 2 = 0, we numerically simu- 
late the cooling and trapping of different nanospheres, 
all with the same initial speed i(0) = 0.10 m/s. The 
time dependent position and the speed of the particles 
are respectively shown in Fig. [2|a) and (b). All parti- 
cles are strongly cooled, and, eventually, trapped in some 
reminiscent small harmonic potential, as can be seen in 
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Fig. [2jc). This takes place within a tenth of a second 
and within a distance of a few millimeter. Both the time 
it takes to be trapped and the distance traveled strongly 
depend on the particle mass, so optomechanical cooling 
in this form allows to sort particles according to their 
mass. The particle speeds shows a sharp jumps around 
v ~ Aj/(2k) = 1.25 x m/s for which we observe 

strongly enhanced cooling. These jumps are obviously 
due to resonances, as can be seen from the amplitude 
in the cooling mode |ox|/|a|, as plotted in Fig. [2jd) . 
The optical amplitude slightly shifts once the particle 
is trapped. This is due to the fact that the particle is 
trapped at a potential minimum for which the optical 
resonance frequency differs slightly from that for a mov- 
ing particle. 

In order to use the mass-dependent cooling and trap- 
ping to weight and sort particle according to their 
mass, we also need to initialize all particles to the same 
speed. This can be achieved by combining optomechan- 
ical damping and heating in our two-mode set up. In 
Fig. [3j we show full numerical results for the parameters 
given_ above, A r = -(4tt x 10 5 + 2k) = -1.4 x 10 6 Hz 
and A 2 = 4-7T x 10 5 — 2k = 1.1 x 10 6 Hz and various ini- 
tial speeds i(0) ranging from 0.02 m/s to 0.16 m/s. The 
results show that all particles reach the same velocity of 
0.1 m/s on the time scale of less than a hundredth of a 
second, which, again, corresponds to a traveled distance 
of few millimeters. 



IV. ANALYTICAL RESULTS 

In order to gain some insight in these features, we ana- 
lytically study the interaction between a single nanopar- 
ticle and the cavity modes. This analysis is general and, 
below, we shall apply it to various parameter regimes and 
applications. 



A. Effective optical force 

In order to derive an expression for the effective op- 
tical force exerted on the particle, we make two sim- 
plifications: (a) if the second term in Eq. pi), which 
describes excitation of the optical mode through the op- 
tomechanical coupling, is much smaller than the first, 
the optical modes can be adiabatically eliminated and 
(b) if a,j << oij, the optomechanical interaction can be 
linearized, i.e. (Aj(z) — Aj)(a>j + ctj) in ^ can be ap- 
proximated by (Aj(z) — Aj)a.j and \a.j + dj\ in (|3]) by 
| ay | 2 + a^a j + a j a* j . 

As a first step to eliminating optical modes, we specify 
initial conditions at t — > — oo, choose Oj(— oo) = and 
formally solve the linearized form of ([2| 



aj(t) 



ICtj 



dt' e (,A '- s)(t -''> A^t/D-A 



(4) 



Substitution of this result in ^ yields a closed descrip- 
tion in terms of z and p, which no longer involves the 
optical modes but, since the integrand in Q depends on 
z(t'), these equations involve integration of z(t') and are 
non-local in time. However, since (Aj(z(t')) — Aj) ~ g, 
it follows from Q that a,j ~ ctjg/(iAj — k), which 
in ([3]) is multiplied by dAj/dz ~ gk so that, to sec- 
ond order in g, if suffices to solve (|3) to zeroth order 
and substitute the result in Q. With |ai| = \a^\ the 
equations of motion (pf to zeroth order describe a free 
particle, so that z(t')~ z(t) + (p(t) /m)(t' — t') and 



Aj(z(t'))-Aj ~ T5 sin ; 



(k (*(*) 



Pit) 



(t'-t) 



where 



the minus sign refers to the first mode. At this level of 
approximation, the modes are not coupled and the op- 
tical forces exerted on the particle remain independent. 
Coupling of the modes only appears in 0(g A ). After sub- 
stitution of Aj(z(t')) in Q, the integral can be evaluated 



ctjgm 

4 \ 2kp — m(Aj + in) 2kp + m(Aj + in) 



a 3=^-h—\7r, ,a 1 - \ — ,x , . J (5) 



and, for |«i| = \a^\ = \ot\, the effective forces are found from ^ 

g 2 \a\ 2 Ajkms'm(2kz) I mKsm(2kz) — (2kp — mAj) cos(2fcz) 



F 



' (2kp — mAj) 2 + m 2 K 2 



;sin(2fcz) — (2kp + mAj) cos(2fcj 
(2kp + mAj) 2 + m 2 K 2 



(6) 



If the particle is not trapped and the momentum change per wavelength traveled is sufficiently small, the effective 
damping force can be obtained by averaging over a wavelength 

F ( = A a \ = 9 2 \a\ 2 km 2 K f 1 1 

i(P) 2-k J Q Z {Z ' P >- 4 \(2kp~mA J ) 2 + m 2 n 2 {2kp + mA 3 ) 2 + m 2 K 2 , 



I 

If, on the other hand, the particle momentum ap- proaches 0, the effective force (|6j) approaches Fj(z) 
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FIG. 3: Velocity as a function of time for particles of size 
r — 50 nm for various initial velocities ranging from vq = 
0.02 m/s to vo — 0.16 m/s in a two-mode setup with one 
damping and one amplifying mode. Particles with velocities 
below v = 0.1 m/s are accelerated while particles with larger 
velocities are decelerated to reach the same stationary velocity 
v — 0.1 cm/s. 



g 2 \a 



2 A ] ksm(4kz)/(2(A 2 + K< 
effective potential 



, which corresponds to an 



g 2 \a\ 2 Aj cos(4kz) 
8(A 2 + k 2 ) 



(8) 



Although the bare optical potentials cancel if |ai| = \ot2\, 
in regime (i), cooled particles are eventually trapped due 
these small additional potentials. Once trapped, the par- 
ticles will perform damped oscillations around their final 
positions. Similarly, in the combined damping and heat- 
ing setup (ii) , the velocity will perform small velocity os- 
cillations, even after it has reached its stationary value. 

The optical forces |7|) and (|6| have sharp resonances 
around p rcs = ±Ajm/(2k). For these values of the mo- 
mentum, the modulation of the intensity of the optical 
field, as perceived by the moving nanoparticle, p/(2mk) is 
resonant with the optical detuning Aj. Such resonances 
give rise to sharp jumps in the velocity of the particle 
and correspond to regimes of strongly enhanced damp- 
ing (or heating). This is illustrated in Fig. [2ja)-(b) and 
in Fig. [3] The width of the resonances is set by the op- 
tical damping rate k according to Ap = Km/ (2k) and it 
follows that for A > n, the sidebands are resolved. For 
p > (p < 0) and a red detuning Aj < 0, the second 
(first) terms in ^ and |7]) are dominant. These terms 
give rise to damping. For blue detunings Aj < 0, it is the 
other way around, and the optical force gives rise to heat- 
ing. Higher-order resonances appear at integer multiples 
of Am /(2k) in higher-orders of g 2 . However, for typical 
experimental parameters, these are highly suppressed. 

Both the average Q and the effective ([6| forces can 
easily be integrated numerically. For the parameters dis- 
cussed above, the results virtually coincide with those 
of our full numerical simulations. This confirms that the 
approximations made are well justified for typical experi- 



mental parameters. As illustrated in Fig. [2^c), somewhat 
larger quantitative differences only appear once the parti- 
cle is trapped by the small reminiscent potentials Q. At 
this point, the assumption that the zeroth order motion 
is that of a free particle is no longer a very accurate one 
and standard-second order optomechanical laser cooling 
theory [T^J [T5] , which assumes harmonic oscillations to 
zeroth order, would give a more accurate effective de- 
scription. In the regime where we are interested in, all 
features are accurately described by the effective force 
The effective damping force ^ almost perfectly de- 
scribes the trend line of the curves. 



B. Mass-dependence 

Having derived the effective optical force, we study 
how the optomechanical damping in regime (i) depends 
on the particle mass. The velocities at which the reso- 
nances occur v les = ±Aj/ (2k) are obviously independent 
of mass, as is also confirmed by the numerical results in 
Fig. [2|a)-(b). The particle mass only appears explic- 
itly in z — p/m and it follows that the optical forces ^ 
only depend on the mass through the momentum, i.e. 
Fj(z,p) = Fj(z 1 mz). However, for actual nanospheres, 
the optomechanical coupling strength g scales with the 
particle volume (see appendix) and, therefore, with the 
mass. Since F oc g 2 , it follows that also Fj oc m 2 so that 
only Fj(z,p)/m 2 is truly independent of mass. Thus we 
find from mz — J2j Fj( z iP) that, for given initial speed 
i(0), the time it takes to be trapped scales with 1/m. 
Since z = J dt z(t), the same is true for the distance trav- 
eled. These conclusions are in obvious agreement with 
the numerical results in Fig. [2Tb). 



C. Optomechanical velocity control 

The combined damping and heating in regime (ii) re- 
quires a little more analysis. In order to gain some in- 
sight, we start from the effective damping force ^ for 
both modes. For the sake of the argument, we assume 
that p > and, for simplicity, we neglect the off-resonant 
terms in Q. In this case, the damping (1) and amplify- 
ing (2) forces are respectively given by 



Gj(p) = T 



g 2 \a\ 2 km 2 n 



4(2kp ± mAj) 2 + m 2 K 2 



(9) 



For equal but opposite detunings A x = — A2, the total 
force Gtot = Gj obviously vanishes. As is illustrated 
in Fig. Ill this is not true for (slightly) different detunings 
Ai = — (A + 5) and A2 = A — S, in which case the total 
force is found as 



6 



G t . 



g 2 \a\ 2 km 2 K 



1 



(2kp — to(A — S)) 2 + m 2 K 2 (2kp — to (A + 5)) 2 + m 2 n 2 



(10) 



This total force 
mA/(2fc) = to(A 2 
is given by 

( dG to 



has a stationary point for p = 
- Ai)/k. The slope around this point 



V dp 



2g 2 \a\ 2 k 2 SK 
m{5 2 + k 2 ) 2 



(11) 



and is negative for 8 > 0, in which case the station- 
ary point is stable. Physically speaking, this means 
that the red-detuned damping laser is further detuned 
than the blue-detuned amplifying laser. Initially too 
fast particles are slowed down, while initially slow ones 
are accelerated to reach the same stationary velocity 
i = pjra = A /(2k), independent of their mass. Assum- 
ing that p < 0, similar arguments lead to the conclusion 
that also p = —mA/ (2k) is a stable stationary point of 
6 > 0. 

From a calculation similar to the one that led to Eq. 
^ , we find a small oscillatory potential for particles with 
velocities around the stationary points 



Vu 



g 2 \a\ 2 8 cos(Akz) 
8(S 2 + k 2 ) 



(12) 



This gives rise to small velocity fluctuations around the 
stationary velocities. For the parameters discussed in 
Sec. Ill, 5 = 2k ~ 1.6 X 10 5 Hz (as for the numeri- 
cal results in Fig. |]), we find g 2 \a\ 2 S/(8(S 2 + n 2 )) ~ 
6.3 x 1CP 24 J, which is almost three orders of magnitude 
smaller than thermal fluctuations at room temperature. 
For a particle with mass to = 1.15 x 10 -18 kg and at a 
speed of 0.1 m/s, it corresponds to velocity fluctuations 
of ~ 5.4 x 10~ 5 m/s. This is in excellent agreement with 
the numerical result shown in the insert in Fig. [3] We 
emphasize that such velocity fluctuations are, at least in 
principle, fully coherent. In this respect they are funda- 
mentally different from thermal fluctuations. 



V. APPLICATIONS 

A. Weighting and sorting 

Combining both settings, we envision the follow- 
ing scheme to sort particles according to their mass. 
Nanoparticles of unknown mass are inserted one by one 
in the optical cavity and confined to the cavity axis. On 
a length scale of a few millimeter, they are either acceler- 
ated or slowed down to a fixed speed of ~ 0.1 m/s using 
the combined optomechanical cooling and heating. Since 
the distance traveled during the velocity rectification de- 
pends both on the initial velocity and on the mass, the 
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FIG. 4: Laser cooling of a 50 nm nanoparticle and i(0) = 
0.25 m/s with (red) and without (black) a compensating po- 
tential. Without the compensating mode, the particle oscil- 
lates rapidly and is quickly trapped so that this type of laser 
cooling becomes much less efficient. Even with a compensat- 
ing potential, the particle is eventually trapped, but the range 
of strong exponential damping is much larger. For parame- 
ters, see main text. 



particles need to be detected when they pass a certain 
transverse plane of the cavity. Once a particle has been 
detected, the detuning of the heating mode is switched 
from the blue to the bare cavity resonance, so that the 
other mode cools the particle until it is trapped, again on 
the length scale of a few millimeter. Both steps easily fit 
within the focal range of a cavity of length L = 0.12 m, 
where the cancellation of the optical potentials is (al- 
most) complete. The distance traveled from the detec- 
tion until the trapping is inversely proportional to the 
mass of a particle and the mass can be determined from 
the position where the particle is trapped. If the effects 
of thermal fluctuations are sufficiently reduced by lower- 
ing the background pressure, the measurement accuracy 
is limited by the position measurements. These are, at 
best, accurate up to A = 10~ 6 m, which amounts to a 
relative error of the order of 10 -3 in the particle mass. 
Switching times are limited by k/2 — 4 x 10 4 Hz, which 
for velocities around 0.1 m/s correspond to position er- 
rors of, at most, a few wavelengths. The effect of velocity 
fluctuations due to the small optical potential ^ amount 
to even smaller relative errors of the order of 10 -4 . 

For the parameters discussed here, mediated interac- 
tions between different particles are weak, and numerical 
simulations on two particles confirm that the small os- 
cillatory optical potential is strong enough the keep the 
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first particle trapped while the second one is inserted 
and slowed down. For beam widths much larger than 
the particle diameters, the cross section of collisions will 
be relatively small and, for just a few particles, all this 
could work within a single cavity. For larger particle 
numbers, collisions become more likely. In the spirit of 
Ref. [TT], once trapped, one could let particles fall down 
into a second storage cavity to avoid collisions. 



B. Transverse cooling towards the quantum regime 

The experimental parameters discussed in Sec. |III| 
were chosen so as to illustrate the resonances, the mass- 
dependent damping and velocity rectification, as shown 
in Fig. [2]and[3j They lead to rather modest overall cool- 
ing. The effective cooling rate is largest if the detunings 
are chosen such that the initial speed of the particle is 
around or just below the resonance and efficient cooling 
only works up to the point where the particle is trapped. 
In this regime of \p\ < Am/ (2k), the effective optical 
damping force Q exerted by a single optical mode on 
a nanoparticle gives rise to an approximately constant 
cooling rate 



7o P t 



2fcV|q| 2 A K 
m(A 2 + k 2 ) 2 



(13) 



For the experimental parameters discussed in Sec. |III[ 
this equals 7 op t — 8.9 x 10 3 Hz and gives rise the the small 
domains with strong exponential cooling Fig. [2^,. In 
principle, the cooling rate can be improved by several or- 
ders of magnitude by using smaller cavities, in which the 
optomechanical coupling is stronger, and higher intra- 
cavity powers. It may be comparable to or even ex- 
ceed the cooling rate in the resolved side band regime 
~ 3x 10 5 Hz, as estimated in pQ. However, larger | ct | 2 <? 2 
also enhances the residual optical potential (jSJ), so that 
one also has to reach a compromise between the range of 
velocities that is cooled and the cooling rate. Smaller cav- 
ities and larger intra-cavity powers are currently operated 
to cool molecules and nanoparticles in quantum interfer- 
ence experiments with nanoparticles and large molecules 
[TU] . In this setup the particles are shot through the cav- 
ity, as illustrated in Figjlja) and the transverse motion 
is strongly cooled during the passage time. 

In order to illustrate the cooling potential of our two- 
mode setup and the improvement compared to standard 
single-mode laser cooling, we consider a smaller cavity 
with L = 0.01 m and a larger intra-cavity power of 
P = 3.2 x 10 2 W. We choose Ai =_-re/\/3 so as 
to optimize the cooling rate ( 13 ) and A = for the 
other, compensating mode. Leaving all other parame- 
ters unchanged, we find k ~ 9.4 x 10 4 Hz, which through 
Ai = —n/sfZ and v rcB — A\/(2k) corresponds to a reso- 
nant speed of 4.3 x 10~ 2 m/s. From we find that the 
particle is trapped for speeds below 8.7 x 10~ 3 m/s, which 
leaves a reasonably large range for exponential cooling. 
In the absence of the compensating mode, the particle is 



trapped in an effective potential of strength g|a| 2 , which 
corresponds to a velocity of 0.24 m/s and leaves no range 
for efficient cooling at all. In Fig|4j we show full numer- 
ical results for both setups. These clearly illustrate the 
enhancement of laser cooling due to the compensating 
potential. 

In principle, the range of velocities can be further 
increased by also compensating the effective potential 
^ or by adiabatically varying the detunings and drive 
strengths. However, eventually, the particle will always 
be trapped in some small optical potential and, in order 
to reach the motional quantum ground state, one has to 
apply resolved side-band cooling as well. The resonances 
and potential cancellation we have discussed in this pa- 
per could be used to enhance pre-cooling in the regime 
where the particle is not yet trapped and the strengths 
may be optimized so as to make the transition between 
the two cooling regimes as fast as possible. 



VI. SUMMARY AND CONCLUSION 

We have studied the classical dynamics of sub- 
wavelength dielectric nanoparticles in an optical cavity 
with two laser-driven modes. The modes where chosen 
such that the optical potentials along the cavity axis can- 
cel and the dynamics is dominated by dissipation. In a 
real system, such cancellation is always an approxima- 
tion, but numerical simulations have confirmed that, at 
least in the regime in which the particle is not trapped, 
a small imbalance of the optical intensities hardly affects 
the dynamics. We have focused on the case in which the 
optical intensities are strong enough to confine the parti- 
cle to the cavity axis, but our results also apply to setups 
in which particles are shot through the cavity orthogonal 
to the cavity axis. In either case, the dynamics is effec- 
tively one dimensional. This setup is somewhat reminis- 
cent of Doppler cooling of atoms but, since nanospheres 
do not have any resolved internal structure, it is also fun- 
damentally different. 

The availability of two resolved optical modes natu- 
rally leads to two different regimes: (i) in which the par- 
ticle motion is either damped or amplified by the optical 
modes and (ii) in which the particle motion is damped 
by one mode and amplified by the other. In case (i), we 
find mass-dependent damping and heating, while regime 
(ii) gives rise to tunable stable stationary velocities. We 
have derived intuitive analytical expressions of the effec- 
tive force on the particle, the effective dissipative force 
and the perturbative optical potential in which cooled 
particles are eventually trapped. The analytical results 
are in excellent agreement with full numerical simula- 
tions. This analytical treatment is accurate up to second 
order in the optomechanical coupling strength g. Since g 
is proportional to the particle mass, it follows that, the 
effective optical force scales with m 2 so that, in case (i) 
and for a given initial velocity, both the time it takes 
for a particle to be cooled and trapped and the distance 
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traveled are inversely proportional to the particle mass. 
This allows for weighting and sorting of nanoparticles 
according to their mass. The measurement accuracy is 
limited by the position measurements which are, at best, 
accurate up to a wavelength. For realistic parameters 
this allows for a relative error on the order of 10" 3 . In 
regime (ii), the combination of heating and cooling leads 
to stationary stable velocities, independent of the parti- 
cle mass. This can be used to initialize different particles 
to the same speed. 

The setup and approach discussed in this paper are 
more generally applicable. Both the mass-dependent 
damping (or heating) and the velocity rectification may 
be optimized for specific needs by adiabatically vary- 
ing the detunings of the laser drives. Thus it should 
be possible to benefit from the resonant cooling over a 
larger range of velocities and to reduce the residual ve- 
locity fluctuations after rectification by another order of 
magnitude. Moreover, we hope that the approach dis- 
cussed here will also prove fruitful in the ongoing efforts 
to create and manipulate coherent quantum states of di- 
electric nanoparticles. A full quantum treatment would 
require a detailed analysis of the effects of thermal de- 
coherence and optical scattering. For the cooling, this 
would be very similar to the analysis in Ref. [T] . Heating 
and, in particular, the combination of damping and heat- 
ing poses additional challenges since it involves relatively 
large speeds and additional fluctuations. At present, it 
is not yet clear how these compete with quantum co- 
herence. Obviously, both heating and combined damp- 
ing and heating may have interesting applications in the 
quantum domain. Heating would, for instance, allow for 
a strong enhancement of quantum fluctuations and the 
combination of heating and cooling could drastically in- 
crease the velocity spread of a quantum wave packet. 
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Appendix A: Optomechanical interaction 

The interaction between a single optical cavity mode 
and a dielectric particle derives from a dependence of the 
optical resonance frequency lo{z) on the particle position 
z along the cavity axis. With the particle inside, the 
displacement field in the cavity is given by D = e E + P, 
where E is the optical electric field and P = eoa p 0(r) 
is the polarization due to the particle. Here a p = 3(e — 
to) /(e + 2eo) is the polarizability of the particle with e the 



permitivity [55] and 0(r) equals 1 for r inside the particle 
and and elsewhere. Under the assumption that \i = 
Ho, so that H = /-i " 1 B, we find from the homogeneous 
Maxwell equations the modified wave equation 



y2 _ (l + a p e (r)) 8 



dt 2 



E = 



(Al) 



with c the speed of light. In the absence of the particle 
(0(r) — > 0), this equation reduces to the standard wave 
equation, which, for given mirror positions and curva- 
tures, can be solved for the complete and orthonormal 
set of mode functions F^ and frequencies ujj — ckj, with 
V 2 Fj = kjFj and where j abbreviates spatial and polar- 
ization mode indices. If the effect of the additional term 



a6(r) in (Al) is sufficiently small, it can be treated per- 
turbatively. To first order, the mode functions are not 
modified. The eigenvalues u> 2 , on the other hand, are 
shifted by — a p J c efar F|(r)0(r)F :) (r), where C denotes 
the volume of the cavity. Since the integrand differs from 
zero only for r inside the particle, the integration can be 
restricted to the volume of the particle. When the parti- 
cle size is much smaller than the optical wavelength, the 
variation of the mode functions Fj(r) over the particle 
volume can be neglected and we find 



WA 1 



aV 

IT 



F,(ro)f 



(A2) 



where r denotes the (center-of-mass) position of the par- 
ticle. For the standing- wave modes Fi^(z) introduced in 
the main text, this reduces to cji^z). 



Appendix B: Full Hamiltonian 



In general, the vector potential corresponding to a 
set of cavity modes Fj is given by A = EjGF( r ), 
from which the canonical momentum density is found as 

II(r,i) = eo A(r,t) = eoEj&WFiW = E^i (<)**(*)• 
From the general Hamiltonian of electrodynamics in the 
radiation gauge i? rad = ^ J ^(I^n - £qC 2 A 1 'V 2 A) 
[27] and the orthogonality of the mode functions Fj, 
we find the cavity Hamiltonian H cav = J^jrf/fao) + 
(eowf (z)/2)£j, which, of course, takes the form of the 
Hamiltonian of a set of uncoupled harmonic oscilla- 
tors. We introduce the normal variables Cbj{z) = 
y/eQUjj{z)/2 £j +i/ yj2tQbj j(z) iTj and their complex con- 
jugates a*(z), so that H cav — J2j u j( z )\ a j( z )\ 2 - We re- 
strict ourselves to two cavity modes and assume that the 
others are sufficiently off-resonant so that they can be 
neglected. As discussed in the main text, the modes are 
separated by a free spectral range of the cavity, which 
we assume is much smaller than both uj^ and oj^ , but 
much larger than any other relevant time scale. It follows 
that the frequency difference {lu^ — c*4 I can be neglected 
with respect to the absolute frequencies lo^ and u>^ but 
that the modes can be addressed independently. The 
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laser drives are modeled by i/drive = Ylj=i i ^j( a j( z ) + 
a*(z))(e l "j + e l ^ ). We move to a frame rotating 
with the drive frequencies a — > a^e i , neglect terms 

rotating with ±2tu^ and assume that g << u> so that 
the z-dependence of a can be neglected, while we keep 



that of the detunings Ai :2 (z) = lo\ . 2 — ^{ ^{z). Thus, 
we obtain the full Hamiltonian in Eq. 0. So far we 
have assumed that z is fixed but this analysis applies to 
z(t), provided that the time dependence is adiabatic with 
respect to the optical time scales. 
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